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Abstract 

The Relativistic Flux Tube model is used to calculate of the Isgur- Wise 
Q>^ . functions describing the exclusive semileptonic decays of B and Bg mesons. 

The light quark mass dependence is investigated and the predicted univer- 
sal function agrees well with the results of lattice simulations, and with the 
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I experimental data. Recent experimental measurements of the B — > D*lvi 

decay distribution yield the CKM element V^c = 0.035 ± 0.001. The IW func- 
. tion slope and second derivative at the zero recoil point are predicted to be 

^ ■ ^'(1) = -0.93 ± 0.04 and f (1) = 1.7 ± 0.1 for a range of light quark massses. 



The importance of including higher derivatives in analyses of experiment is 
emphasized. 



1 Introduction 



In the heavy quark hmit (mass of the heavy quark mq oo) all the non-perturbative, 
strong interaction physics for semileptonic B Dlvi and B D*Wi decays (and 
corresponding B^ decays), can be parametrized in terms of a single universal func- 
tion |T|], known as the Isgur-Wise (IW) function. Knowledge of the IW function 
is essential in all calculations that lead to numerical values of branching ratios of 
semileptonic decays and of the Vcb element of the CKM matrix. A reliable calcu- 
lation of the IW function from first principles is still beyond our power, although 
recent efforts using lattice QCD are very encouraging. A variety of different meson 
models make specific predictions for the IW function. Since these model results 
vary widely, the IW prediction can be considered to be a sensitive gauge of model 
validity. 

The relativistic flux tube (RFT) model shows promise to provide a realistic 
description of all meson states. The RFT model is in essence a description of dy- 
namical confinement [0-0. For slowly moving quarks rigorous QCD relativistic 
corrections clearly demonstrate that the scalar confinement potential picture 

is incorrect 0, ^. On the other hand, the RFT model dynamics is consistent with 
both spin-dependent |TD|, ^ and spin-independent |^ QCD expectations. In 



this paper we examine the RFT model predictions for IW functions. We find here 
that the spinless light degrees of freedom model is in excellent agreement with both 
experiment and with fundamental QCD expectations. 

We begin by a brief review of the RFT formalism. In Section 3 we outline the 
general calculation of the IW function and its derivatives at the zero recoil point. 
Our numerical results and conclusions are given in Section 4. Here we evaluate the 
RFT predictions for the IW function and examine its dependence on the light quark 
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mass. We determine the Vet CKM element and branching fractions for various 
exclusive semi-leptonic decays and test the shape of the predicted IW function. 
Finally we note that the RFT prediction for the slope at the zero recoil point is in 
good agreement with present lattice simulation results. 



2 RFT model in the heavy quark hmit 

In the limit ttiq oo the quantized equations of the spinless RFT model with light 
quark mass nig are P, || 



/(/ + !) 1 

— = ^{Wr,j±v^} + a{rJ{v^)} , (1) 

1 TTT. T (i r arcsinv_L, , 

H, = H-mQ = {Wrn±} + 7:{r, -} + Vir). (2) 

z / v± 

Here Hg is Hamiltonian of the light degrees of freedom (LDF), with the short range 
potential 



K 



and 



V{r) = -- , (3) 



Wr = ^Pl + ml , (4) 

7x = -j=L^, (5) 

1 arcsinf_L 1 , . 

f{v±) = — ( ) . (6) 

Av± v± 7_L 

It is assumed that v± is Hermitian, and the symetrization (curly brackets in (|I|) and 
(0)) then yields a Hermitian Hamiltonian. 

The solution of the energy eigenvalue equation 

Hg<^ = Eg<!> (7) 
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required the development of methods to ehminate the a-priori unknown operator v± 
and also to address the solution of operator equations in which functions of both pr 
and r appear. Our approach H, || employs the Galerkin method ||12| in which 
a truncated basis set is used to transform the operator equations into finite 

matrix equations. The resulting matrix equations describe the motion of a light 
quark with radial momentum and perpendicular velocity v± relative to a static 
heavy antiquark. In order to solve for the eigenvalues of the Hamiltonian, one first 
solves the angular momentum equation (|I|) for the v± matrix, and then diagonalizes 
the Hg matrix (0) as discussed in detail in references P, |[. 

Since there is no spin-orbit coupling in the above model, we can write our wave 
function (describing the LDF in the meson rest frame) as a product of the orbital 
wave function and the light quark spinor as 

,(0) _ ^(0) (^\A0) ^-iEnt 



*i^i(x) = $L%,,(x)0^i^e--^^ (8) 

Here, Hq represents all quantum numbers of the LDF, and the superscript denotes 
the rest frame of a meson. 



3 Decays Hi Hcli^i and the IW functions in the 
spinless RFT model 

For heavy-light mesons, the spin of the heavy quark sq and the spin of the light 
antiquark decouple as and in the limit uiq — > oo they are conserved separately 
by the strong interactions. Because of this, hadrons containing a single heavy quark 
can be simultaneously assigned the quantum numbers {sQ,ms^) and {sg,ms^), and 
matrix element of the hadron current describing semileptonic decay can be expressed 
by the free heavy quark current and the IW function. 



Except for a trivial kinematical factor, the IW function $,{uj) is defined |l4 



as the overlap of the wave functions describing the light degrees of freedom in the 
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two mesons, i.e. 

e(c^) = {"^mi^H,) ■ (9) 

Here, if the meson velocity in the lab frame is v, the wave function describing its 
LDF is 

^H,ix') = S{^r)¥|ll{x) , (10) 

with x' = A~^(v)x being the lab frame, x the rest frame of the particle, and S{\) is 
the wave function Lorentz boost. 

Since the IW function is Lorentz invariant, we can choose any frame to calculate 



it. Particularly convinient is the modified Breit frame |13| , |T4|, where the two parti- 
cles move along the z-axis with equal and opposite velocities. Denoting velocity of 
He as V, by the use of ( p!OD the overlap takes the form 

= {<I>H,\<PH,) I d'x'^U^')^H-Sx%'=0 (11) 

= (02il'5^(v)^(-v)|0g)) / rfV<I>gf (x+)5t(v)5(-v)$g)(x„)|,=o . 

In this expression x+ and X- denote the rest frames of He (moving in the +z direc- 
tion) and Hi (moving in the -z direction), respectively. Using the fact that Lorentz 
boosts are real, i.e. S^'iy) = S(y) = 5'~^(— v), the boost factors cancel out. The 
overlap of spin part of the wave functions will just give us S h^, which we sup- 
press in the following, and we are left with 

i^H-J^H,) = J d^x'^gf (X+)$S)(X_)|.=0 . (12) 

Denoting v = |v| and 7 = :^7j=f , and using 

x±\t'=o = A(±v)a;'|t/=o = (TIVz' , x' ,y' ,'jz') , (13) 

we obtain 

(vl/^jvl/^j = ld'x'^^;^^Ji~^vz',x',y',^z')^P^i+W,x',y',^z') 

= J d^x'^^I^Jix', y\ iMty. y', T^Oe-'^^^^''^' , (14) 
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since the energies of the hght degrees of freedom are the same in both hadrons. 
Finally, after rescaling the z' coordinate (2;' ^-z')? renaming integration variables, 
and noting kinematical identities 



w + 


1 


2 




1 


u; + 


1 



7=,/^. (15) 

(16) 



valid in the Breit frame, the IW function becomes [14 



e(^) = ^ / rf^x<|.gt(x)<|.g)(x)e-2^^^^^ . (17) 

By spherical symmetry, and since we are interested here only in s-waves, the 
orbital wave function has the form 

$aOo(x) =i?„o(r)^ . (18) 
v47r 



For the calculation of the IW function we need to evaluate integral of the form 

I c?3xi?:,o(r)i?,o(r)e-2^^'^^ . (19) 



Using 



00 



e 



-ikz 



in 



Y^(2l + l){-tyUkr)^I^^Yio, (20) 



1=0 

and orthonormality of the spherical harmonics, one easily obtains the IW function 
in the form 

«.) = -l^(,„(2£,/^.)), (21) 

where 

{A) = / dr r^R^,o{r)A{r)R^o . (22) 
Jo 
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The result (|21| ) can also be used to obtain derivatives of the IW function at the 
zero recoil point [uo = 1). Explicit results for the first two such derivatives are 

\ + \El{r')). (23) 

(24) 



^(1) 



Since all integrals in these expressions are positive, we recover the limit ||l4| 

1 



and also obtain 



r(l) < 

r(i) > 



2 ' 
1 



(25) 
(26) 



Similar limits can be found for the third and higher order derivatives. 



4 Results 

The parameters appearing in the spinless RFT model are the string tension a, the 
short range potential constant k, and the quark masses. In our calculation, we fix 
the tension from the universal Regge slope a' ~= 0.8 GeV^"^, 

a = ~ 0.2GeV^ . (27) 

This is consistent with the value found from analyses of heavy onia spectroscopies 
[1T2| . We also fix the light quark mass rriu^d, as the quality of our fit depends only 
weakly on its value. The remaining parameters k, rris, nic and are varied to best 
account for the 6 spin averaged states [§, |^. As an example, in Table |1| we show 
the results of our fit when rriu^d was fixed at 0.3 GeV. Once the parameters of the 
model are determined, the wave functions and LDF energies Eg are known and the 
IW function can be calculated with the aid of 

The IW function, as computed from ([21|), is displayed in Fig. |l|. The solid and 
dashed curves correspond to rriu^d = and 300 MeV respectively. We observe that 
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with increasing light quark mass the two IW functions differ by only a few percent 
at uj = 1.5. 

One can expand the IW function about zero recoil point as discussed at the end 
of the preceeding section. Using ( ]23|) and (p^) we have explicitly evaluated the first 



and second derivatives over the usual range of light quark masses rriu^d as shown in 
Figs. 1^ and 0. These derivatives vary by less than ten percent over the range of 
niu^ and we conclude that 

^'(1) = -0.93 ± 0.04 , (28) 
^(1) = 1.7 ±0.1. (29) 

In Fig. ^ we compare the power series approximation about the zero recoil point 
a; = 1 to the full IW function. We observe that the power series expansion is not 
very convergent and that a significant error is encountered even including cubic 
terms in — 1 near the high end of the u range. In particular, if one truncates the 
series to include only the linear term, then in fitting to the data the slope will be 
overestimated, i. e. the true slope will be more negative. 

The explicit values of the derivarives will provide an excellent test of the RFT 
model through lattice simulation calculations. It has been proposed [jl5| that a direct 



evaluation of higher derivatives of the IW function at zero recoil point is possible by 
numerical QCD methods. 

In Fig. ^ we show the ratio of the IW functions for and B semileptonic 
decays. Using light quark mass values rus = 513 MeV and niu^d = 300 MeV, from 
the fit of Table |l|, we see that the IW function for the two decays differ by less than 
by 1.5%. We also note that ^Bsi^) — Cb{^)- 

A comparison of our IW function ^b{'-^) with the recent ARGUS |TB| and CLEO 



II data fT^ (shown in Fig. P) yields the Vcb element of the CKM matrix 

Kb = 0.035 ±0.001 (30) 
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(the of this fit was 0.88 per degree of freedom). The RFT prediction for the 
uj dependence of the IW function is in excellent agreement with the data. Finally, 
using expressions given in flSl for differential widths of B and decays, we calculate 



corresponding branching ratios to be 



Br{B - 


DW) 


= 1.79 


Br{B - 


> D*liy) 


= 5.01 


BriBs - 


. DJu) 


= 1.83 


Br{B, - 


> d:w) 


= 5.11 



1.53ps 

'_[B_, 

1.53ps' 
1.53ps' 



70 . 



(31) 
(32) 
(33) 
(34) 



Agreement with the available data [|T9] for B D (1.8 ± 0.4)% and for B ^ D* 
(4.5 ± 0.4)% is satisfactory. We also find that 



Br{B DW) ^ Br{B, D^Ip) , 



and 



BriB D*W) ^ BriB. 



d:w) 



(35) 



(36) 



as expected from the heavy-quark symmetry. 

Since the IW function is normalized at the zero recoil point to ^(1) = 1, much of 
the model dependence is distilled into the slope at zero recoil ^'(1). Various models 
have yielded a wide range of slopes from — | > .^'(1) > —2. Our result for the 
slope in the RFT model is from (E3) 



4(1) = -0.93 ±0.04, 



(37) 



where the error is estimated from the slope dependence on the light quark mass 

is 



shown in Fig. The recent CLEO II experimental slope value [|17 

4(1) = -0.82 ±0.12 ±0.12 , 



(38) 
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which is consistent with our predicted value. This value was obtained by assuming 
(^(cj) ~ 1 + ^'(l)(a; — 1). As we have discussed the slope obtained by this truncation 
will be overestimated, i. e. the actual value should be more negative. 

The RFT model must of course agree with experimental data, but it also must 
be consistent with fundamental results from QCD. The lattice simulation method 



provides several recent results for ^'(1). The calculation of Bernard et al. ||21| gives 

^^(1) = -1.24 ±0.26 ±0.33 , (39) 



while that of Mandula and Ogilvie 1 22|| yields 

4(1) ^ -0.95 



and the UKQCD collaboration result is 



-1 2+°-'^ 

^•^-0.3 



(40) 



(41) 



The agreement among these results is a gratifying demonstration of the consistency 
of the RFT model with QCD. 
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TABLES 



Table 1: Heavy-light spin averaged states. Spin-averaged masses are calculated in 
the usual way, by taking | of the triplet and | of the singlet mass. 



state 


spectroscopic label spin-averaged 


theory 


error 






2^+1 L J mass (MeV) 


(MeV) 


(MeV) 


cu, cd quarks 












D (1867) 
D* (2010) 


0" 
1- 




IS (1974) 


1975 


1 


Di (2424) 


1+ 




IP (2424) 


2425 


1 


cs quarks 












Ds (1969) 
D* (2110) 


0- 
1- 


'So' 
'S, 


1 c (or\7r:\ 
ij yz,u 1 ) 


2074 


-1 


Ds, (2537) 


1+ 




IP (2537) 


2536 


-1 


bu, hd quarks 












B (5279) 
B* (5325) 


0- 
1- 


'So' 
'S,_ 


IS (5312) 


5311 


-1 


bs quarks 












Bs (5368) 


0- 


'So' 


IS (5409) 


5410 


1 






mu,d 


= 300 MeV (fixed) 










rus 


= 513 MeV 










rric 


= 1292 Me\^ 










nib 


= 4628 MeV 










a 


= 0.2 GeV^ (fixed) 










K 


= 0.516 
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FIGURES 



Figure 1: IW function for B D^*Hvi decays for mu,d = and mu,d = 0.3GeV 
(dashed line). 

Figure 2: Slope of the IW function at zero recoil for B decays as a function of mu,d- 

Figure 3: Second derivative of the IW function at zero recoil for B decays as a 
function of mu^^- 

Figure 4: Power series expansion about oj = 1 (zero recoil point) truncated at various 
number of terms. The solid curve is the full IW function for B decays. In this figure 
we have taken mu,d = 0.3 MeV. 

Figure 5: Ratio of and for the particular choice of rriu^d = SOOMeV. The 
mass of the strange quark was rus = 513MeV as a result of the fit of Table 0. We 
observe the two IW functions differ by less than two percent over the kinematically 
allowed range of uj. 

Figure 6: Comparison of the IW function ^si'^) calculated from the RFT model for 
the particular choice of rriu^d = 300 MeV with ARGUS [|16] and CLEO II data [jl 



for the decay B^ —>■ D*^l ui. For the sake of clarity, error bars are shown only for 
the CLEO II data. 
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Figure 1 
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Figure 6 
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